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The  high  cost  of  modern  weapons  systems  force  increasingly 
stringent  requirements  on  structural  integrity  and  durability 
of  new  military  hardware.  Prediction  of  durability  in  storage 
and  service  environments  requires  therefore  an  increase  in 
sophistication  which  necessitates  a strong  interaction  between 
experiment  and  analysis. 

This  report  is  Part  III  of  an  investigation  into  the 
moisture  effects  in  carbon  fiber  reinforced  epoxy  composites. 
Specifically,  it  deals  with  the  prediction  of  moisture  sorption 
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natural  environments.  Such  knowledge  is  essential  for  the 
analysis  and  prediction  of  mechanical  property  changes. 
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INTRODUCTION 

A commitment  has  been  made  by  the  Navy  to  use  carbon  fiber 
reinforced  composites  in  future  Naval  aircraft  as  structural 
materials.  Stringent  requirements  on  structural  integrity  and 
durability  prior  to  aircraft  production  provide  further  impetus 
to  fuller  understanding  and  resolution  of  the  moisture/thermal 
environmental  effects  and  their  prediction  in  a real  service 
environment. 

This  report  is  Part  III  of  an  investigation  into  the  effects 
that  moisture  has  on  carbon  fiber  reinforced  epoxy  composites. 

In  Part  I and  Part  II  we  discussed  diffusion,  and  changes  in  the 
matrix  dominated  strength  properties  as  a function  of  the  combined 
effect  of  moisture  concentration  and  temperature.  This  report 
discusses  the  prediction  of  moisture  uptake  as  a function  of  time 
and  the  corresponding  internal  distribution  in  carbon  fiber  compos- 
ites in  a real  outdoor  or  service  environment.  Such  knowledge 
is  essential  for  the  analysis  of  the  concomitant  changes  in  the 
elastic  properties  which  will  be  discussed  in  a forthcoming  report, 
and  we  also  believe  that  it  is  essential  for  a theoretical  explanation 
and  prediction  of  the  reversible  strength  degradation. 

EXPERIMENTAL 

For  the  experimental  measurements  of  the  moisture  equilibrium 
concentrations  in  the  composite  matrix  as  a function  of  relative 
humidity  and  for  the  measurements  of  composite  diffusion  coeffi- 
cients as  a function  of  temperature,  we  refer  to  Part  I of  this 
report  m. 


DISCUSSION 

The  subject  to  be  discussed  here  is  the  diffusion,  uptake, 
and  internal  distribution  of  moisture  in  a composite  exposed  to 
real  outdoor  environments  under  service  and  storage  conditions. 

Other  degradative  mechanisms  such  as  UV  radiation  damage  or  thermal 
spiking  and  mechanical  fatigue  will  not  be  considered.  There  is 
no  doubt  that  these  latter  degradation  mechanisms  may  be  superimposed 
in  real  service  environments,  and  cannot  be  neglected;  however,  these 
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mechanisms  depend  more  on  the  structural  application  and  have  to 
be  ascertained  independently.  We  are  presently  investigating  the 
effect  that  thermal  and  mechanical  cycling  has  on  the  diffusion 
of  moisture  into  the  composite.  These  results  shall  be  reported 
at  a later  date. 

Also,  the  analysis  of  the  effect  that  moisture  (and  its  internal 
distribution)  has  on  composite  elastic  constants  will  be  discussed 
in  a forthcoming  report.  The  knowledge  of  the  internal  moisture 
distribution,  however,  is  essential  for  such  an  analysis. 

A.  Prediction  of  Moisture  Absorption  and  Internal  Moisture 
Distribution  In  Organic  Materials. 

If  one  would  know  exactly  the  future  environmental  history 
of  a given  composite  element,  one  could  predict  with  a high  degree 
of  certainty  what  the  moisture  level  and  internal  distribution  of 
moisture  would  be  at  any  time.  All  that  is  required  is  to  solve 
Fick's  second  equation  (1) 


Where  c = concentration,  t = time,  and  Dt,  D2,  D3  are  the  principal 
diffusion  ficients,  the  magnitude  and  direction  of  which  are 

describr  rotational  ellipsoid.  As  we  have  stated  in  Part 

I of  th  1 the  diffusion  coefficients  (Di,D2,D3)  depend 

on  the  volume  fraction  (Vf) , on  the  temperature,  and  perhaps 

on  the  cure  concentration  (depending  on  the  matrix). 

For  a composite  plate,  where  the  edge  effects  can  be  ignored, 
equation  (1)  reduces  to  the  one  dimensional  equation  (2). 


where  Dz  is  the  diffusion  coefficient  vertical  to  the  composite 
plane. 


The  boundary  conditions  are  assumed  such  that  moisture  penetrates 
from  both  sides  of  the  sheet,  i.e.,  the  concentration  on  the  surface 
of  the  composite  is  equal  to  the  solubility  (=  equilibrium  concen- 
tration) of  moisture  in  the  composite,  and  the  initial  interior 
moisture  distribution  is  known. 


In  a real  outdoor  environment  the  boundary  conditions  change 
continuously  because  the  relative  humidity  (RH)  changes.  And,  since 
the  equilibrium  concentration  changes  with  the  partial  pressure 
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of  moisture  in  the  air  so  does  the  surface  concentration  on  the 
sheet.  Over  a daily  or  seasonal  RH  cycle  there  will  be  absorption 
and  desorption,  but,  the  long-term  net  effect  will  be  an  absorption 
till  equilibrium  has  been  reached  with  the  surroundings.  Therefore, 
we  can  say  that  the  final  moisture  equilibrium  depends  on  some  kind 
of  RH  average. 

Also,  the  temperature  changes  continuously  and  so  does  the 
corresponding  diffusion  coefficient,  which  may  be  described  by  an 
Arrhenius  equation 


D = Dq  exp  (-E/RT ) 


(3) 


where  D0  is  a constant,  E is  the  activation  energy  of  diffusion, 

R is  the  ideal  gas  constant  and  T is  the  absolute  temperature  in  °K. 

Assuming  that  we  have  experimentally  determined  the  moisture 
equilibrium  concentration  as  a function  of  RH  and  temperature  (as  in 
Figure  1)  and  the  diffusion  coefficient  (as  a function  of  temperature 
as  in  Figure  2,  and  if  necessary,  as  a function  of  concentration) , 
and,  that  we  know  what  the  environment  will  be  (or  has  been) , we  are 
now  in  a position  to  calculate  the  moisture  absorption  and  internal 
moisture  distribution  as  a function  of  time. 

The  problem  can  be  solved  in  two  ways: 

a)  by  some  averaging  of  temperature  and  humidity  conditions 
or 

b)  by  a finite  difference  approach  where  the  environmental 
temperature  and  humidity  changes  are  introduced  in  discrete 
finite  time  intervals  and  by  solving  the  diffusion  equation 
(which  has  been  converted  to  a finite  difference  equation) . 

Both  procedures  will  be  described  in  some  detail  in  the 
following. 

B.  Solution  of  the  Diffusion  Equation  by  a Temperature  and 
Humidity  Averaging  Procedure! 

a)  Series  Solution. 

If  the  RH  and  the  temperature  in  the  environment  had  the  known 
constant  values  = M(RH]_)  , and  the  initial  moisture  distribution 
had  the  uniform  value  Co»  then  the  moisture  uptake  and  internal 
distribution  would  be  given  immediately  by  equations  (4a  and  4b)  £2} . 


{,2  j 3!  Crank , "Mathematics  of  Diffusion"  Oxford  (1956)  . 
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jp  » 1 - f ^ — 2 exp  [--D(2n+l)2Ti2t/4h2] 


(4a) 


n=0 


(2n+l) 


epe;  ■ 1 - ? J UHTlf  exP  t -D(2n+l)2,2t/4h23  cos  ,4b) 


n=o 


where  Mt  = amount  of  moisture  absorbed  after  a time  tf  M°°  amount 
of  moisture  absorbed  after  an  infinite  time,  i.e.,  equilibrium 
concentration  or  solubility;  h = half  of  the  thickness  of  the 
composite  sheet,  and  D corresponding  to  T]_. 


Figure  3 shows  a reproduction  of  the  solution  from  reference  T2], 
in  dimensionless  coordinate.  X/h  = 0 corresponds  to  the  composite 
midplane,  and  (C-Cg)/(C.  - C~)  = 1 corresponds  to  the  constant  surface 
concentration  in  equilibr iumwith  the  surrounding  RH.  The 
(Dt/h2)l/2  parameter  gives  a master  curve  (Figure  4)  which  holds 
for  all  materials  following  Fickean  diffusion.  If  Mt  and  M°°  is 
replaced  by  the  percent  moisture  concentration  after  time  t and 
by  the  equilibrium  percentage  concentration  respectively,  one  obtains 
directly  the  moisture  uptake  (in  %)  as  a function  of  time  and  the 
problem  is  solved. 


If  we  knew  the  "right"  temperature  and  humidity  averages  of 
an  environment  where  temperature  and  humidity  changes  continuously, 
we  would  also  have  the  average  solution  for  such  an  environment. 


In  other  words  the  "right"  averages  Tkav  and  RHkav  are  those 
that  would  give  the  "same"  results  had  the  composite  been  exposed 
to  constant  weather  conditions  with  temperature=Tkav  and  humidity  = 

RHkav. 


From  inspection  of  equation  (3)  we  see  that  one  cannot  use 
the  arithmetic  average  temperature  Tav  since  the  diffusion  coeffi- 
cient is  exponentially  dependent  on  the  temperature,  i.e.,  for  a 
given  temperature  cycle  the  diffusion  rates  at  the  higher  temperatures, 
T > Tav,  are  more  dominant  than  those  at  the  lower  temperatures 
T < Tav»  Therefore,  one  has  to  use  an  average  diffusion  coefficient 
which  is  obtained  by  converting  the  temperature  curve  to  a diffusion 
coefficient  curve  of  which  the  average  integral  value  (Dkav)  is 
used  for  the  calculation.  (Dkav  will  be  called  the  kinetic  average 
diffusion  coefficient  and  the  corresponding  temperature  Tkav  will 
be  called  the  kinetic  average  temperature.)  From  equation  (3)  it  is 
also  easily  seen  that,  the  larger  the  temperature  fluctuation  is, 
the  more  is  the  deviation  of  Tkav  from  Tav,  and  Tkav  is  always  larger 
or  equal  to  Tav 
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In  a similar  way  we  can  show  that  it  is  not  permissible  to 
use  the  arithmetic  average  of  the  relative  humidity  cycle  since 
the  RH  values  over  a daily  cycle  are  generally  lowest  when  the 
temperatures  are  highest,  i.e.,  in  a sample  with  a uniform  internal 
distribution  corresponding  to  the  arithmetic  average  humidity  there 
would  be  a higher  rate  of  desorption  during  the  hot  part  of  the 
cycle  than  absorption  during  the  cold  part,  thus  there  would  be 
a net  desorption.  On  the  other  hand,  the  rate  of  diffusion  depends 
also  on  the  concentration  gradient,  which  is  highest  during  the 
cold  part  of  the  cycle.  This  effect  may  be  even  more  pronounced 
if  the  equilibrium  concentration  increases  non  linearly  at  high 
relative  humidities  (>85%  RH) . The  determination  of  RHfcav  will  be 
discussed  below  under  b. 

Thus,  one  can  establish  correction  factors  Ft  and  Frh  for  the 
temperature  and  relative  humidity,  which,  multiplied  with  the 
arithmetic  average  temperature  and  humidity  gives  the  kinetic 
average  temperature  and  humidity  respectively.  The  higher  the 
temperature  fluctuation  over  a daily  or  seasonal  cycle  the  larger 
are  Ft  and  F^.  Ft  and  Frh  will  therefore  be  material  constants  for 
a specific  climatic  (geographic)  area  (or,  more  generally,  of  a 
military  hardware  service  cycle  in  a given  geographic  area) . 

These  factors  (Ft  and  Frh)  can  be  easily  determined  for  various 
geographic  locations  or  for  climatic  categories  from  meteorological 
data  and  now  it  becomes  very  easy,  for  engineering  purposes,  to 
estimate  moisture  absorption  in  organic  materials  without  the  use 
of  high  speed  computers.  (An  example  of  such  an  engineering  problem 
is  described  in  reference  £3J) 

b.  Determination  of  Fm  and  F_,„ 

x Kn 


Since  Tav  and  RHav  can  be  readily  obtained  from  almanacs  or 
encyclopedias,  the  knowledge  of  Ft  and  Frh  for  various  climatic  cate- 
gories may  be  all  that  is  required  for  estimates  in  engineering 
diffusion  problems  for  outdoor  conditions. 

11 

For  the  composite  material  under  consideration  it  is  assumed 
that  for  any  temperature  value  T one  can  obtain  the  corresponding 
diffusion  coefficient  for  moisture  D = D(T)  (e.g.,  from  equation  (3)), 
and,  for  each  relative  humidity  value  H,  one  can  obtain  the  corre- 
sponding solubility  (equilibrium  percent  by  weight  moisture)  M = M(H) 

(e.g.,  by  interpolation  of  tabulated  experimental  values).  We  present 
an  averaging  procedure  for  obtaining  values  for  Tfcav  and  RH)<av  given 
8 evenly  spaced  temperature  and  relative  humidity  readings  over 
a 24-hour  period.  The  values  of  T|<av  and  RHkav  will,  of  course, 
depend  on  the  material  being  considered;  specifically  they  will  depend 


ur  J.  M.  Augl  "Prediction  of  Moisture  Absorption  in  the  5"54 
HIFRAG  Discarding  Rotating  Band  for  Various  Climatic  Environments." 
NSWC/WOL/TR  76-167 
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on  the  functional  relationships  D (T)  and  M(H) . From  the  description 
of  this  case  it  should  be  clear  how  to  apply  this  procedure  to  average 
any  tabulated  weather  data  (such  as  data  at  1-hour  or  3-hour  intervals 
over  a month  or  a year) . 

Let  the  8 temperature  and  humidity  values  be  denoted  by  , 

T2»  . ..,  Tq  and  by  H^,  H2,...,Hgr  pespectively.  Let  D be  a fixed 
positive  number  (for  example,  set_D  = D (293^15  °K) ) and  let  T be 
the  temperature  corresponding  to  D (so  D(T)=D).  Now  the  change 
in  moisture  concentration  throughout  a composite  plate  after  3 hours 
of  being  at  temperature  T and  whose  surface  during  this  time  is 
maintained  at  moisture  concentration  M(H)  is  exactly  equal  the  change 
in  moisture  concentration  throughout  the  same  plate  after  3 D(T)/D 
hours  of  being  at  the  temperature  T and  whose  surface  during  this 
time  is  maintained  at  moisture  concentration  M(H) . This  follows 
quickly  from  changing  the  time_scale  in  the  diffusion  equation 
Mt  = DMxx  (since  if  x = D (T) t/D  then  MT  = DMXx) • 

Thus  one  can  approximate  the  effect  of  the  weather  data  given 
by  the  eight  3-hourly  readings  Ti,  Hi  (i=lfi. .,8)  by  the  following: 
let  the  composite  be  exposed  for  xi  = 3D(T]_)/D  hours  at  temperature 
T with  surface  moisture  concentration  MJHi) , followed  by  exposure 
for  X2  = 3D(T2)/D  hours  at  temperature  T with  surface  moisture  concen- 
tration M(H2) , and  so_on,  ending  with  exposure  for  x8  = 3D(Tq)/D 
hours  at  temperature  T with  surface  moisture  concentration  M(Hs) . 

We  approximate  the  latter  situation  by  the  case  where  the  composite 
is  exposed  for  x = (x^  + t2  + • . . + x8)  hours  at  temperature  T 
and  with  surface  moisture  concentration 


1 


i'i 


Mkav  - {T1  M(H1}  + T 2 m(H2)+...+t8  M(H8))/x. 


(5) 


As  explained  above,  exposure  for  x hours  at  temperature  T and  with 
surface  moisture  concentration  M^av  is  exactly  equivalent  to  exposure 
for  24  hours  with  diffusion  coefficient 


D,  = D 
kav 


x/24 


(6) 


and  with  surface  moisture  concentration  Mkav*  Let  Tkav  be  the 
temperature  corresponding  to  Dkav  (so  Dkav  = D(Tkav))  and  let  RHkav  be 
the  relative  humidity  corresponding  to  Mkav  (so  Mkav  = M (RHw av) ) . 

Then  exposure  of  the  composite  material  to  the  weather  data  *Ti,  h^} 
may  be  approximated  by  the  situation  of  exposure  to  the  constant 
weather  condition  temperature  = Tkav»  relative  humidity  = RHkav 
for  the  same  time  period  (which  in  this  example  was  24  hours) . 
Numerical  results  obtained  using  this  procedure  are  given  in 
Tables  1 and  2. 
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c.  Solution  of  the  Diffusion  Equation  by  a Finite 
Difference  Method. 


From  meteorological  stations  one  can  obtain  weather  data 
(including  temperature,  RH,  wind  and  sky  conditions)  in  three-hourly 
intervals  over  many  years  on  magnetic  tapes.  These  data  are  quite 
suitable  to  be  read  directly  by  a high  speed  computer  and  can  be 
used  to  solve  the  diffusion  equation  by  a finite  difference  method 
in  discrete  three-hourly  time  steps  over  a period  of  several  years, 
and  assuming  that  such  a multi-year  weathercycle  remains  the  same 
in  the  average  if  one  takes  multiples  of  it,  the  process  can  be 
repeated  as  often  as  necessary  to  obtain  the  final  answer. 

The  finite  difference  method  can,  as  well,  be  used  to  solve 
the  diffusion  equation  when  monthly  or  seasonal  or  yearly  kinetic 
average  temperatures  and  relative  humidities  have  been  obtained. 

The  boundary  data  and  diffusion  coefficient  in  the  moisture  diffusion 
equation  are  then  changed  only  every  month  (or  season  or  year) 
instead  of  every  three  hours,  and  hence  the  numerical  solution 
procedure  takes  much  less  computer  time.  In  the  case  of  Narmco 
5208/T300  epoxy  composites  it  will  be  seen  that  the  results  of 
3-hourly  data  can  be  excellently  approximated  using  monthly  averages 
even  with  a relatively  thin  composite  thickness  (Figure  6) . 

We  now  describe  the  finite  difference  method  used  to  solve 
the  diffusion  equation  for  moisture  in  a plate  of  composite  material 
with  thickness  2h.  The  diffusion  equation  is 


Mfc  = (D(t)Mx)x  for  0<x<2h,  t>0 

where  M(x,t)  denotes  the  percent  by  weight  of  moisture  in  the 
composite  (at  location  x and  at  time  t) , and  D is  the  diffusion 
coefficient.  Since  D depends  on  temperature,  and  the  variation 
of  the  temperature  with  time  is  assumed  to  be  given,  one  may  consider 
D to  be  directly  a function  of  time.  Both  surfaces  of  the  composite 
plate  are  assumed  to  be  exposed  to  the  same  relative  humidity  condi- 
tions, so  the  moisture  concentration  M(t,x)  in  the  plate  will  be 
symmetric  about  the  midplane  x=h,  and  it  is  thus  sufficient  to  solve 
the  following  diffusion  problem: 


= ( D ( t ) Mx ) for  0<x<^h,  t>0 

M(x,t=0)  = M° ( x ) for  O^x^h 

M (x=0 , t)  = g(t)  and  Mx  (x  = h,t)  = 0 for  t>0 

13 


(7a) 

(7b) 

(7c) 
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where  the  subscripts  t and  x symbolize  partial  differentiation  with 
respect  to  time  and  direction  vertical  to  the  plane  of  the  plate,  and 
where  M°(x)  is  the  given  initial  moisture  concentration  within  the 
plate  (generally  M°=0) , and  g(t)  is  the  known  surface  moisture 
concentration  which  is  equal  to  the  solubility  of  moisture  in  the 
composite  corresponding  to  the  (given)  relative  humidity  at  time 
t.  The  values  D(t)  and  g(t)  are  determined  from  the  given  weather 
data  (or  from  Tkav  and  RH|<av  values  which  were  obtained  from  the 
weather  data) . 


The  diffusion  problem  (7)  can  be  solved  numerically,  using 
the  standard  implicit  finite  difference  method  which  is  described 
below.  Let  I be  a positive  integer,  define  xq=0  and  xj=h,  and  pick 
points  xi,...,xi~l  such  that  xQ<xi< . . . <xi~1<xi . For  i=l,...,I 
define  A^  = x^  - xj-1,  and  choose  a time  step  At.  Let  N denote 
the  number  of  time  steps  to  be  taken  (the  final  time  at  which  an 
approximation  to  the  moisture  concentration  in  the  composite  is 
obtained  is  then  N*At) . The  numerical  solution  value  approximating 
M(xi,n«At)  will  be  denoted  by  M.  (for  i=0,...,I  and  n=0,...,N). 
Define  M.  = M°(Xi)  for  i=0,...,I,  and  define  m"  = g (n*At)  for  n=0, 
. . . , N . if  values  M?,  i=0,...,I,  have  been  obtained  approximating 
the  moisture  concentration  inn£ljie  plate  at  time  tn  e n*At,  then  the 
approximate  solution  values  M?  approximating  the  moisture  profile 
at  the  next  time  level  tn+-*-  = (n+l)At  can  be  obtained  by  solving 
the  following  linear  system  of  I equations: 


M n+1  n 
Mi  Mi 

At 


2D(tn+I) 

VAi+l 


,n+l 


/ 1 1 \ n 

— ‘I1!  ^~i) Mi 


for  i = 1, . . . , 1-1 


n+1 


Mn+1 

Mi+1 

Ai+1 


(8a) 


and 


m n+1  Mn 
Mi  ~ Mi  = 2Di 
At 


r— 1 (<;  - Mf1) 


.n+1 


at  i = I 


.n+1 


(8b) 


Note  that  (8b)  results  from  defining  M,  + , = mt_i  an<^  AI  + 1 - A j -:md 
letting  i=I  in  (8a) . The  average  percent  by  weight  moisture  in 
the  composite  at  a given  time  t**  is  defined  to  be 


i f M(x,tn)dx 


(91 
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and  was  numerically  approximated  by 


t 5 


(10) 


Further  details  on  the  numerical  method,  including  choices  of  the 
spatial  mesh  ( A i ) and  the  time  step  (A^)  and  the  method  of  solution 
of  the  linear  equations  (8)  are  given  m Appendix  A. 

D.  Eight  Categories  of  Climatic  Conditions 

Guidelines  for  realistic  considerations  of  climatic  conditions 
in  research,  development,  test  and  evaluation  have  been  published 
by  the  Army  under  Army  Regulation  AR  70-38,  July  1969.  This  report 
divided  the  earth's  climatic  conditions  into  eight  climatic 
categories : 

Category  1,  wet-warm 
Category  2,  wet-hot 

Category  3,  humid-hot  coastal  desert 

Category  4,  hot-dry 

Category  5,  intermediate  hot-dry 

Category  6,  intermediate  cold 

Category  7,  cold 

Category  8,  extreme  cold 

A reproduction  of  the  world  map  indicating  these  climatic  areas 
from  AR  70-38  is  shown  in  Figure  5.  For  all  categories,  a represen- 
tative daily  temperature  and  humidity  cycle  is  given  for  "operational" 
and  "storage  and  transit"  conditions.  We  have  calculated  the  D|<av 
and  RH|^av  of  Narmco  5208/T300  for  all  eight  categories  for  the 
operational  and  storage  conditions  (Table  1)  and  compared  these 
values  with  the  Dav  values  corresponding  to  the  average  temperature 
(Tav)  and  the  average  relative  humidity  RHav.  The  ratios  of  these 
values  give  the  Fr  and  Frh  factors.  The  way  these  correction  factors 
should  be  used  is  to  look  up  the  yearly  temperature  and  humidity 
averages  of  a particular  environment,  then  determine  to  which 
climatic  category  it  will  belong.  Then  the  averages  are  multiplied 
by  the  corresponding  Ft  and  FRR  factors.  The  Tkav  and  RH|<av  values 
thus  obtained  are  used  as  pseudo  constant  temperature  and  humidity 
conditions . 

It  should  be  remembered  that  the  F and  FRH  factors  depend  on  the 
climatic  condition  and  on  the  material  (i.e.,  on  the  slope  of  the 
log  D vs.  1/T  curve  and  on  the  equilibrium  concentration  curve) . 

The  effect  of  the  climatic  category  on  these  factors,  however,  is 
much  more  pronounced,  as  can  be  seen  from  Table  2 where  we  compared 
the  Tkav,  RHkav  and  the  Ft  and  FRh  of  Narmco  5208/T300  CF  composite 
with  Nylon  612  which  is  a totally  different  material. 
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The  Tuav  can  be  considered  equal  for  all  practical  purposes, 
but  also  the  RH^av's  are  surprisingly  close.  It  appears  therefore 
that  these  Ft  and  Frh  correction  factors  can  be  used  for  moisture  dif- 
fusion problems  in  most  organic  materials  as  a first  approximation. 

In  tropical  areas,  such  as  the  Panama  Canal  Zone,  there  is  little 
difference  between  seasons,  and  the  temperature  fluctuations  between 
day  and  night  are  so  small  that  T|<av  at  Tav  and  RH^y  a*  RHav 

E.  Limitation  of  the  Averaging  Procedure 

a)  Moisture  Absorption. 

An  averaging  procedure  may,  by  its  very  nature,  be  insensitive 
to  fluctuations.  The  fluctuations  in  moisture  content,  caused  by 
seasonal  temperature  and  ; midity  changes  in  the  environment,  are 
obviously  stronger  in  tb  omposite  plates  than  in  thick  ones. 

While  a calculation  usin^  uriable  boundary  conditions  at  short 
time  intervals  will  give  x alistic  fluctuations  in  moisture  uptake, 
an  averaging  procedure  will  give  only  smooth  curves. 

Examples  are  given  in  Figures  6 and  7.  Figure  6 shows  a computer 
printout  of  moisture  uptake  for  a 6-ply  Narmco  5208/T300  composite 
laminate  if  it  had  been  exposed  at  Washington,  D.C.,  under  climatic 
condition  between  1972  to  1975  (this  3 yearly  cycle  was  repeated). 

The  solid  line  was  obtained  by  supplying  the  computer  with 
3 hourly  temperature  and  humidity  data  and  the  stars  represent  an 
averaging  over  each  month  of  the  year. 

One  can  see  that  an  averaging  over  monthly  periods  is  almost 
as  good  as  3-hourly  data  inputs.  It  requires  however  considerably 
less  computer  time. 

The  repeat  cycle  of  36  months  is  clearly  visible  which  reflects 
that  there  are  visible  differences  also  in  the  yearly  climatic 
changes. 

The  thicker  the  composite,  the  less  pronounced  are  these 
fluctuations  as  shown  in  Figure  7.  These  curves  show  the  moisture 
uptake  of  6,  12,  18,  24,  30,  36,  42  and  48  ply  laminates  (.01397  cm 
per  ply) . 

The  smooth  curves  were  obtained  by  using  a 3-yearly  average 
while  the  oscillating  curves  were  obtained  from  monthly  averages. 

For  24-ply  laminates  it  is  quite  sufficient  to  use  long-term 
average  data  of  T|<av  and  RH|<av. 


Also,  as  one  might  expect  in  temperate  zones  the  largest  amount 
of  moisture  is  absorbed  during  the  summer  months,  where  the  average 
temperatures  are  highest.  Figures  8 and  9 show  computer  printout 
for  6-  and  18-ply  laminates  respectively  if  such  composites  would 
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have  been  exposed  starting  with  January,  May  and  September  1972. 

Such  details  can,  of  course,  not  be  observed  if  yearly  averages 
would  have  been  used.  These  differences  are  significant  only  during 
the  initial  moisture  uptake. 

b)  Internal  Moisture  Distribution 

Kinetic  average  temperature  and  moisture  conditions  can  be 
used  to  obtain  the  time  history  of  the  internal  moisture  distri- 
bution within  a composite  exposed  to  variable  weather  conditions. 

The  results  of  using  monthly  kinetic  averages  (see  Table  4)  to  calcu- 
late the  moisture  distribution  within  6 plies  (half  thickness 
= .04191  cm.)  of  Narmco  5208/T300  with  fiber  volume  fraction  = .7 
agree  very  well  with  the  values  obtained  by  running  the  corresponding 
3 hourly  weather  data,  except  (as  to  be  expected)  for  a thin  layer 
near  the  surface  which  reflects  what  the  last  several  hourly  readings 
happen  to  have  been:  see  Figure  10  in  which  the  solid  lines  join  the 
values  obtained  at  the  mesh  points  of  the  finite  difference  method 
used  to  solve  the  diffusion  problem  with  month  by  month  kinetic 
averages,  and  in  which  the  values  obtained  when  using  3 hourly 
weather  data  are  plotted  as  individual  points. 

For  this  thin  laminate,  the  result  of  using  the  constant 
moisture  boundary  condition  and  diffusion  coefficient  corresponding 
to  the  3 year  kinetic  average  of  the  weather  data  (Mkav  = *82975, 

Dkav  = 5.556E-11)  is  not  a good  approximation  to  the  actual  moisture 
uptake  (see  Figure  11  in  which  the  lines  are  as  in  Figure  10  while 
the  individual  points  were  obtained  using  the  3 year  kinetic  average) . 
In  both  Figure  10  and  11  a solid  line  is  plotted  at  .82975  = Mkav» 
and  curves  are  plotted  at  28  and  45  months.  When  using  the  monthly 
kinetic  averages,  the  weather  readings  for  1972-1974  were  continually 
repeated,  and  in  this  thin  composite  the  moisture  profiles  became 
virtually  3 year-periodic.  After  28  months,  the  maximum  and  mini- 
mum moisture  profiles  (when  using  month  by  month  kinetic  averages) 
occur  at  28  and  45  months,  respectively,  and  all  other  profiles 
lie  in  between  these  two  curves. 

For  a thicker  laminate,  using  the  3 year  kinetic  average 
produces  a reasonably  good  approximation  (Figure  12).  Since  the 
kinetic  average  was  taken  over  a 3 year  period,  one  would  expect 
the  best  approximation  to  occur  at  3 year  intervals  (i.e.  at  multi- 
ples of  the  period  over  which  the  averaging  was  done) . Here  the 
approximation  (away  from  the  boundary)  at  other  times  is  almost 
as  good.  In  Figure  12  moisture  profiles  after  1,2,3,7,8,9,15,16, 

17  years  in  36  plies  (half  thickness  = .25146  cm.)  of  Narmco  5208/ 
T300  with  fiber  volume  fraction  = .7  are  given.  The  lines  were 
obtained  using  month  by  month  kinetic  averages  while  the  points 
were  obtained  using  the  3 year  kinetic  average. 
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TABLE  4 

MONTH  BY  MONTH  KINETIC  AVERAGE  MOISTURE  SOLUBILITY  AND  DIFFUSION 
COEFFICIENT  FROM  WEATHER  DATA  AT  WASHINGTON,  D.C. 

FOR  RACH  MONTH  THE  3-H0URL7  '-'FATHER  OATA  HAS  BEFN  CONVERTED  (USING  THE 
MFTHOD  OFSCRIRFO  IN  THIS  RFpohT)  TO  FOOT V ALF.N T AVERAGES  FOR  FXPOSU°F 
OVER  A TO  DAY  T I ME  PERIOD  (MAP^CO  S2n8/TT00  •«'ITh  FIBER  VOLU^F  FRACTION 
EQUAL  .7).  THUS  EXP^SURF  Of  THIS  COMPOSITE  TO  THE  wFAThER  CONDITIONS 
OUR  I MG  EACH  MONTH  TS  AOPPOX  IMATFL.Y  EOUI  VALENT  TO  EXPOSURE  TO  THF 
CORRFSPONDING  CONSTANT  AWFRaGF  CONDITIONS  FOR  A DURATION  OF  TO  HAYS. 
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F.  Additional  Corrections  for  Exposure  to  Direct 
Sun  Radiation. 

So  far  we  have  taken  temperatures  reported  in  weather  data. 

This  would  mean  that  the  composite  is  completely  protected  from 
sun  (or  even  sky  radiation) . 

Unpainted  carbon  fiber  composites  that  are  exposed  to  direct 
sun  radiation  can  become  considerably  warmer  than  the  surrounding 
air  temperature,  and  therefore,  the  diffusion  rate  is  also  higher 
than  one  would  calculate  from  the  temperatures  of  the  surrounding 
air.  It  is  important  to  estimate  the  error  in  predicting  the 
moisture  absorption  in  a composite  if  one  only  uses  the  temperature 
or  humidity  data  from  a particular  outdoor  exposure  site. 

For  this  purpose  we  have  exposed  four  CF-composite  panels 
outdoors. 

Each  panel  was  mounted  horizontally  on  a piece  of  polypropylene 
foam  and  had  a thermocouple  in  its  center.  Three  panels  were  exposed 
to  direct  sun  radiation  (#1  uncoated,  #2  painted  white,  #3  bonded 
to  a thin  aluminum  foil) . One  reference  sample  #4  was  shielded  from 
the  sun  radiation  by  placing  it  inside  a wooden  box,  open  at  the 
opposite  side  of  the  sun.  The  temperatures  were  recorded  simul- 
taneously (day  and  night)  for  several  weeks. 

Four  typical  temperature  profiles  are  shown  in  Figures  13  to 
16.  These  represent  a clear  sunny  day;  a clear  sunny  but  windy  day; 
a hot,  humid,  hazy  day;  and  a cloudy,  rainy  day. 

The  observed  differences  are  quite  remarkable.  The  differences 
in  temperature  between  the  unprotected  samples  covered  from  and 
exposed  to  sun  radiation  can  be  40-50°F.  The  white  painted  sample 
was  somewhat  cooler  than  the  sample  coated  on  both  sides  with  an 
aluminum  foil.  Even  the  unexposed  samples  absorbed  enough  scattered 
radiation  to  become  slightly  warmer  than  the  reference  thermocouple 
not  connected  to  a composite  (see  Figure  13-16) . A clear  day  with 
blue  sky  showed  a greater  temperature  difference  than  a hazy  day. 

A clear  but  windy  day  had  a somewhat  cooling  effect.  And  the  least 
difference  was  observed  on  a cloudy,  rainy  day. 

If  one  had  calculated  a kinetic  average  diffusion  coefficient 
Dkav  (corresponding  to  the  air  temperature) , then  the  calculated 
amount  for  the  moisture  absorption  of  the  sun  exposed  samples  could 
be  off  by  a factor  of  2.  Thus  an  airplane  with  CF  composites  inside 
a hanger  or  outside  on  a runway  would  absorb  different  amounts  of 
moisture.  For  refined  estimates  such  correction  factor  should  not 
be  neglected  in  predicting  real  environmental  effects. 

Thus,  if  unpainted  carbon  fiber  composites  are  exposed  to  direct 
sun  radiation,  then  one  also  should  know  the  sky  conditions  and 
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the  fractional  times  the  composite  is  exposed  to  these  conditions. 
The  diffusion  coefficient  corrected  for  sun  radiation  is  given  by 
(ID 


Dkavr  " Dkav  <2*05fcs  + 1*64fcw  + ^^hh 


+ 1# 14fcr  + V 


(ID 


where  fcs,  fcw,  fhh»  fCr'  and  fP  are  the  fractional  times  the 
composite  is  exposed  to  clear  sunny  sky,  clear  sunny  sky  but  with 
wind,  humid  and  hazy  conditions,  cloudy  rainy  conditions,  and  fully 
protected  from  sun  and  sky  radiation.  The  constant  2.05,  1.64,  etc., 
were  taken  from  the  experimentally  determined  values  for  unpainted 
panels.  (It  should  be  remembered  that  in  the  Dkav  values  the 
correction  factor  for  the  climatic  category  is  already  included.) 

So  that  the  whole  procedure  may  be  shown  schematically  by 


look  up  ! 

l 

RH 

av 

| corr . for_  1 

RH  * frh 
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^kav 

1 corr.  for_  j 
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| climatic  | 
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| radiation  ^ 
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For  storage  conditions  where  only  the  air  temperature  governs  the 


materials  temperature,  this  reduces  to  D 


kavr 


= D 


kav* 


If  the  temperature  fluctuations  over  daily  and  seasonal  periods 
are  small,  then  Dkav  = Dav  and  RHkav  = RHav. 


G.  Effect  of  Coatings  as  a Moisture  Barrier. 


While  organic  coatings  give  a definite  protection  against  UV 
radiation,  we  do  not  expect  that  they  will  be  of  great  benefit  as 
moisture  barriers.  Some  of  the  highly  fluorinated  or  chlorinated 
hydrocarbons  such  as  Teflon,  polyvinylidene  fluoride,  and  polyvinyli- 
dene  chloride  may  perhaps  reduce  the  rate  of  moisture  permeation 
into  the  composite  slightly,  but  will  not  prevent  it. 

In  order  to  estimate  the  effect  of  a particular  coating  one 
needs  to  know  the  diffusion  coefficient  or  the  permeabilities  of 
the  coating  film.  The  moisture  permeability  (P)  is  given  by  (12) 


P = SD 


(12) 


where  S = equilibrium  concentration  or  solubility  of  moisture,  and 
D = diffusion  coefficient. 

The  effective  change  in  the  flux  of  moisture  through  a unit 
surface  between  an  uncoated  and  a coated  sheet  of  material  can  be 
easily  visualized  by  considering  a steady  state  flux  through  such 
a sheet  where  on  one  side  the  humidity  is  kept  at  some  constant 
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value  and  on  the  other  side  it  is  kept  essentially  at  zero.  Let 
us  call  substance  (1)  as  the  coating  and  substance  (2)  as  the  pro- 
tected substrate.  At  the  surface  of  the  coating  the  moisture 
concentration  is  Cq  = Si  (where  Si  is  the  equilibrium  moisture 
concentration  or  solubility) . The  concentrations  on  both  sides  of 
the  interface  (i)  are  not  the  same,  but  are  given  by  the  distri- 
bution coefficient  C^/c}  = k which  have  the  same  ratio  as  the 
solubilities,  i.e.,  S2/Si  = k.  The  internal  concentration 
distribution  of  moisture  is  then  indicated  by  the  sketch  below. 


S2 

q2  ~D2  L2 

.u  ■ ^ 

«12  * -2 

q12  = ^2 

q2  S2 


where  q2  is  the  flux  without  coating  and  q12  Is  the  flux  throughout 
the  coated  plate;  Cj,  and  C~  are  the  concentrations  at  the  inter- 
face, Si  and  S2  are  the  solubilities,  Li  and  L2  are  the  corresponding 
thicknesses  and  Di  and  D2  are  the  diffusion  coefficient  of  substances 
(1)  and  (2)  respectively.  Since  at  the  steady  state  the  flux  is 
the  same  throughout  the  slab,  one  obtains 


Using  permeabilities  P-^  = S^  and  P2  = S2  D2,  one  obtains  further 


P1S1 

S1L1 


siLi 


P2C2 

S2L2 


. Pj  k Sj  Lj  S2  L2 

2 L1(P1S2L2  + P2kSlLl^ 
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Now  the  reduced  flux  due  to  the  coating  can  be  calculated  from  the 
ratio. 


q12  _ ^1  P1  k S1  L2 

q2  S2  (P1  S2  L2  + P2  k S1  Ll) 


P2  L1 
1 + — — 

P1  L2 


or 


kD„  L. 
1 + — - — 


(13) 


Table  3 shows  values  of  the  reduced  flux  through  the  slab  for  various 
ratios  of  permeabilities  and  ratios  of  coating  to  substrate  thick- 
nesses. It  is  quite  obvious  that  a coating  with  a very  low 
permeability  (about  1/1000  of  that  of  the  substrate)  would  be  required 
to  give  a useful  protection,  or  the  thickness  of  the  coating  has 
to  be  increased  considerably. 


The  only  possible  candidates  that  might  be  considered  to  have 
some  effect  are  polymers  such  as  teflon,  polyvinylidene  chloride, 
polyvinylidene  fluoride  and,  perhaps,  a few  highly  fluorinated 
polymer  films.  Of  course,  an  organic  coating  will  still  have  a 
protective  effect  against  ultraviolet  radiation. 


Metal  coating,  on  the  other  hand,  would  be  sufficient  even 
in  very  thin  films  to  protect  the  composite  from  moisture.  In 
addition,  metal  coatings  would  also  give  added  lightning  strike 
protection. 


H.  Engineering  Examples. 

Example  1. 

Determine  the  moisture  uptake  (in  percent)  of  Narmco  5208/T300 
epoxy  composites  with  varying  numbers  of  plies. 

The  thickness  for  a commercial  ply  of  5208/T300  in  a fabricated 
laminate  is  approximately  .01347  cm  (.0055  inches).  The  diffusion 
coefficients  of  a composite  with  a fiber  volume  fraction  of  .70 
has  been  measured  or  can  be  calculated  1 , and  the  temperature 
dependence  is  given  in  Figure  2.  The  solubility  profile  of  moisture 
in  Narmco  5208  is  given  in  Figure  1.  Using  these  data  and  the 
3-hourly  weather  data  of  Washington,  D.C.  (from  1972-1974)  the 
following  three  yearly  averages  were  determined  as  described  under  C. 

Tay  = 58 . 19°F  Tkav  = 62 . 31°F 

°Tav  = 4,818  E-11  cm2/sec  Dkav  = 5.556  E-ll  cm2/sec 
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Table  3.  Relative  change  of  the  rate  of  diffusion  of 
moisture  through  (or  into)  a composite  plate 
(qi2/q2)  as  a function  of  relative  coating 

thickness  and  permeability  ratios. 


Lj/L 


2 


/pl  \ 

0.001 

0.005 

0.01 

0.05 

0.1 

1000 

.50 

.167 

.091 

.02 

.01 

500 

.67 

.285 

.167 

.038 

.02 

100 

.51 

.67 

.50 

.167 

.091 

50 

.95 

.80 

.67 

.285 

.167 

10 

.99 

.95 

.91 

.67 

.50 

5 

.995 

.975 

.95 

.80 

.67 
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RH  = 64.15% 

av 

MRHav  = *8259% 

Ft  = 1.0707 

Where  MRHav  and  MRHkav  are  the  equilibrium  concentrations  of  moisture 
in  the  composite  corresponding  to  64.15%  and  64.45%  of  relative 
humidity,  the  other  symbols  are  the  same  as  defined  before. 

Figure  7 shows  a computer  printout  of  the  percent  moisture 
uptake  of  6,  12,  18,  24,  30,  36,  42  and  48  ply  Narmco  5208/T300 
composite.  The  smooth  lines  were  obtained  by  the  averaging  method 
using  3-yearly  Tkav  and  RHkav  values  and  the  irregularly  shaped  lined 
were  obtained  by  using  monthly  kinetic  averages  in  the  finite 
difference  method. 

Example  2. 

Prepare  a diagram  that  gives  the  fraction  of  the  equilibrium 
concentration  of  moisture  in  a laminate  as  a function  of  number 
of  plies  (or  thickness)  and  time  for  the  same  composite  as  described 
in  Example  1.  The  same  experimental  data  were  used  as  in  Example  1. 

The  results  is  shown  in  Figure  17  which  allows  a quick  estimate 
for  the  moisture  uptake  of  a Narmco  5208/T300  (Vf  = .7)  composite  for 
the  climatic  conditions  of  Washington  D.C.  for  any  number  of  plies. 

If  the  diffusion  coefficient  can  be  considered  as  concentration 
independent,  then  this  graph  is  applicable  for  any  average  relative 
humidity  (RHkav)  (with  the  understanding  that  the  equilibrium  concen- 
tration is  RH  dependent) . For  instance,  a 20-ply  5208/T300  CF 
composite  would  take  6,  16,  39,  76  and  118  months  to  reach  25,  40, 

60,  80  and  90  percent  respectively  of  its  equilibrium  moisture 
concentration  (which,  for  Washington,  D.C.  is  .8275%  by  weight). 

Such  engineering  diagrams  can  be  readily  given  for  any  desired 
environment. 

I . Expansion  to  other  Organic  Materials  and  to  Military 
Hardware  in  Service  Environment. 

The  approach  described  in  this  report  can,  of  course,  be  easily 
expanded  to  other  organic  materials.  It  is  necessary  to  measure 
the  moisture  equilibrium  concentration  profile  and  the  diffusion 
coefficient  as  a function  of  temperature.  More  complex  shapes  of 
components  can  be  usually  broken  down  into  simpler  configurations 
or  the  general  diffusion  equation  in  two  or  three  dimensions  has 
to  be  solved  by  finite  difference  methods. 

For  service  environments  a realistic  mission  profile  is  required. 
Let  us  qualitatively  discuss  an  example. 


RHkav  - 64.45% 


M 


RHkav 


.8297% 


F^  = 1.0046 
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What  is  the  moisture  absorption  of  an  aircraft  wing  skin  made 
of  CF  composite?  In  addition  to  the  skin  thickness,  the  composite 
diffusion  coefficient  as  a function  of  temperature,  the  moisture 
equilibrium  concentration  profile  and  the  local  Tfcav  and  RHk^v  one 
requires  knowledge  about  a projected  aircraft  history.  For  instance, 
fraction  of  time  standing  unprotected  on  the  airfield,  fraction 
of  time  in  the  hangar,  fraction  of  time  in  the  air,  fraction  of 
time  under  supersonic  skin  heating. 

High  altitude  subsonic  flight  will  contribute  very  little  to 
diffusion  because  of  the  low  temperatures.  Supersonic  flights  may 
contribute  significantly  if  a.)  the  composite  is  thin  and  b.)  the  skin 
temperature  is  high  (see  exponential  effect  of  T on  D) . However, 
it  is  expected  that  the  major  contribution  comes  from  the  fraction 
of  time  the  aircraft  is  standing  unprotected  in  the  field. 

The  program  described  in  this  report  can  be  easily  modified 
to  incorporate  such  conditions. 

CONCLUSIONS 


1.  Two  analytic  methods  have  been  described  to  model  the 
moisture  uptake  and  internal  distribution  in  fiber  reinforced  organic 
composites  in  real  outdoor  environments.  Required  for  the  calcula- 
tions are  the  diffusion  coefficients  of  the  composites  as  a function 
of  temperature  and  the  moisture  equilibrium  concentrations  as  a 
function  of  the  relative  humidity. 

2.  The  first  method  is  simple  and  can  be  applied  without  the 
use  of  a high  speed  computer.  It  uses  average  temperatures  and 
average  humidities  combined  with  correction  factors  that  are  specific 
for  certain  climatic  categories  or  geographic  areas  and,  to  a lesser 
degree,  dependent  on  the  material.  This  method  gives  smooth  curves 
that  average  out  seasonal  changes. 

3.  The  second  method  is  exact  but  requires  one  to  three  hourly 
weather  data  over  a period  of  several  years  and  it  also  requires 
access  to  a high  speed  computer,  it  assumes  that  the  data  (of 
several  years)  will  be  representative  for  the  prediction. 

4.  Further  corrections  were  introduced  if  the  temperature 
effect  of  direct  sun  radiation  is  of  interest. 

5.  The  program  can  be  readily  expanded  to  include  various 
profiles  of  military  hardware  mission  environments. 
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SYMBOLS 

CQ,  Clf  C[.  various  expressions  of  moisture  concentration 
D diffusion  coefficient 

Dav  diffusion  coefficient  corresponding  to  the  average  temperature 

D.  diffusion  coefficient  corresponding  to  the  kinetic  average 

av  temperature 

D,  diffusion  coefficient  corresponding  to  the  kinetic  average 

avr  temperature  including  sun  radiation  effects 

Fpj,  factor,  which  multiplied  with  the  average  relative  humidity, 
gives  the  kinetic  average  humidity 

Ft  factor,  which  multiplied  with  the  average  temperature,  gives 
the  kinetic  average  temperature 

f fractional  time  of  exposure  to  cloudy-rainy  conditions 

f fractional  time  of  exposure  to  clear,  sunny  sky  conditions 

C b 

f fractional  time  of  exposure  to  clear,  sunny  sky,  with  windy 

cw  conditions 

f fractional  time  of  exposure  protected  from  direct  sun 

" radiation. 

h half  of  the  composite  plate  thickness 

k distribution  coefficient 

1,L  thickness  of  composite 

M moisture  concentration 

moisture  concentration  at  time  t 

equilibrium  moisture  concentration 

M.  moisture  equilibrium  concentration  corresponding  to  the 

a kinetic  average  relative  humidity. 
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SYMBOLS  Cont. 

P permeability 

q flux 

RHav  average  relative  humidity 

RHkav  kinetic  average  relative  humidity 

S solubility 

t,  t time 

T temperature 

Vf  volume  fraction  of  fibers  in  the  composite 
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EQUILIBRIUM  CONCENTRATION  OF  MOISTURE  IN  NARMCO  5208  RESIN  AS  A FUNCTION  OF 
TEMPERATURE  AND  RELATIVE  HUMIDITY 
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FIG.  3 CONCENTRATION  DISTRIBUTION  (MASTER  PLOT)  AT  VARIOUS  TIMES 
IN  THE  SHEET  -•»<  X < h WITH  INITIAL  UNIFORM  CONCENTRATION  C 
AND  SURFACE  CONCENTRATION  Cv  NUMBERS  ON  CURVES  ARE  VALVES 
OF  Dt/h2  . 


f 
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-Onn  8-000  16.000  24-000  32-000  40-000  48-000  56-000  64-000 

MONTHS 

FIG.  8 PERCENT  MOISTURE  UPTAKE  IN  A 6 PLY  NATMCO  5208/T300  CARBON  FIBER  COMPOSITE  CALCULATED  FOR 

WASHINGTON  D.C.  CLIMATIC  CONDITIONS  (MONTHLY  AVERAGES).  STARTED  AT  VARIOUS  TIMES  OF  THE  YEAR. 
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INTERNAL  MOISTURE  DISTRIBUTION  IN  A 6 PLY  5208/T300  COMPOSITE  (Vf  = 0.7)  AFTER  VARIOUS  EXPOSURE  TIMES, 
COMPUTED  BY  USING  MONTHLY  AVERAGES  (SOLID  LINES)  AND  3-HOURLY  READINGS  (POINTS) 


KINETIC  AVERAGE 
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FIG.  11  INTERNAL  MOISTURE  DISTRIBUTION  IN  A 6 PLY  5208/T300  COMPOSITE  (Vf  = 0.7)  AFTER  VARIOUS  EXPOSURE  TIMES, 
COMPUTED  BY  USING  MONTHLY  AVERAGES  (SOLID  LINES)  AND  3 YEAR  AVERAGE  (POINTS) 
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FIG.  16  TEMPERATURE  VARIATION  OF  CF  COMPOSITES  ON  A CLOUDY  AND  RAINY  DAY.  (AUG.  27,  19761 
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SOLVE  DIFFUSION  EQUATION  ON  INTERVAL  0 TO  XF INAL=  2 . OOOOOOOE-O 1 
NUMBER  OF  X- 1 NTFRV ALS=  10  TOTAL  NUMBER  OF  TIME  STEPS*  4800 

Tl^t  STEP  DT=  1.29b0O00E*05 

OXELE=  • I000000E-02  .3000000E-02  .3000000E-02 

.3000000E-02  .1000000E-01  .1000000E-01  .2000000E-0I 

•5000000E-01  .5000000E-01  .5000000E-01 

T=  o.  0=  1 .OOOOOOOE^OO  JCYCL=  0 

END  OF  MONTH  NUMBER  0 AVE  HDY  CONDITION  DURING  MONTH*  0. 

X=  0.00000  .00100  .00400 

UH=  0.  0.  0. 
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FIG.  19  SAMPLE  COMPUTER  OUTPUT 
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Appendix  A. 

Numerical  Solution  of  the  Moisture  Diffusion  Equation 
a.  Finite  Difference  Equations. 

We  first  give  a derivation  for  the  finite  difference 
equations  (8)  used  to  obtain  an  approximate  solution  to  the  diffusion 
equation  (7).  Define  tn+1  = (n+l)At.  The  following  approximations 
are  used: 


(Al) 

(A2a) 


(A2b) 

•i] 

(A2c) 


Equation  (Al)  accounts  for  the  left  side  of  (8a)  while  (A2c)  with 
(A2a)  and  (A2b)  yield  the  right  side  of  (8a).  Because  of  the  boundary 
condition  Mx(h,t)=0,  or  more  basically,  because  the  solution  of 
the  original  problem  is  symmetric  about  the  midplane  x=h,  one  may 
define  the  approximate  solution  value  Mj  + 1 at  the  point  XI+i  (with 


A-l 


•1 
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XI+1  = Xx  + A,  = h + A ) to  be  Miti  = mt-i*  The  choice  of  X . 
gives  A-  , = A..,  and  tnen  equation  T8a)  Isaefined  for  i=l  ana 
reduces1?©  (8b1. 

b.  The  Linear  Equations. 

The  structure  ofn|he  linear  system  of  I equations  given  by 
(8)  for  determining  M.  (i=l,...,I)  will  now  be  described.  Recall 
that  M?  (i=0,...,I)  and  M~  (n=0,...,N)  are  defined  by  the  given 
initial  data  M^(x)  and  byuthe  given  boundary  data  g(t),  respectively. 

Multiplying  both  sides  of  (8a)  by  (Ai  + Ai+1)/2  produces 


'i  + 4i+l  / t.n+l  ,.n\  D/,n+i\  k-J  (l  , 1 \ n+1  , Mm 

2fit  7 ' 1 ) FT  M‘  *771. 


( (A3) 


Equation  (A3)  at  i=l  has  the  form 


. Mn+1  . . „n+l  _ „ 

A^  f ^ M^  + A^  2 M2  F ^ 


where 


A1+A2 

1,1  2A  t 


♦ D<tn+1>(i7  * 4)'  Al-: 


. -pit"*1). 


(A4a) 


and  F 


(A1+A2)M1  D(tn+1)  n+i 

- '5TZ + — a M_ 


1 2A  t 


”0 


Equation  (A3)  for  i in  the  range  2 through  1-1  has  the  form 


**<*-*  nT- 1 + Ai,i  M"+1  + Ai.i+i  * Fi 


where 


A = 

i,i-l  Aj 


1/1 


Ai  + Ai+1 
2A  t 


i , i+1 


-p(tntl) 


+ D(tn+‘)(^  + fe) 


4+1 


. (Ai  + Ai+l)..n 
d Fi  " 2A t Mi . 


( A4b) 
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Multiplying  both  sides  of  (8b)  by  Aj/2  produces  the  equation 


where 


. un+l  , A „n+l 

^1  I— 1 ^1  — 1 ^1  I ^1 


= F , 


1,1-1 


-D(t 


n+] 


A M 

and  FI  = "lAtT 


I, I 
n 


h_  + D(tn+1) 

2A  t Aj 


( A4c) 


Considering  A to  be  an  I by  I matrix  (entries  of  A which  have  not 
been  specified  in  equation  (A4)  are  defined  to  be  zero) , and  consid- 
ering Mn+1  and  F to  be  vectors  of  length  I,  the  equations  (A4)  , 
which  are  equivalent  to  (8) , have  the  form 


AM 


n+1 


F. 


(A5) 


The  matrix  A and  vector  F in  (A5)  are  those  which  result  from 
solving  (7)  by  using  the  finite  element  method  (fully  implicit  in 
time)  with  linear  elements  and  with  lumping  of  the  mass  matrix  (by 
using  the  trapezoidal  rule)  DQ,  and  this  is  the  linear  system  of 
equations  which  was  solved  in  our  computer  program  to  obtain  Mn+1. 

c.  Evaluation  of  A and  F in  Equation  A5. 

The  entries  of  A and  F depend  on  the  time  step  At,  the  spatial 

mesh  (i.e.,  the  A^)  , the  value  D(tn+^),  the  value  Mq  + \ and  the  known 

approximate  solution  values  Mn  at  the  previous  time  level.  Choices 
for  At  and  the  Aj  will  be  discussed  later.  The  diffusion  coefficient 

D tn+l  depends  on  the  temperature  at  time  tn+^  while  m!?+‘'‘  is  the 
equilibrium  concentration  (solubility)  of  moisture  in  the  composite 
corresponding  to  the  relative  humidity  at  tn+l.  Thus  D(tn+1) 
and  Mq  depend  on  the  weather  data  and  on  the  particular  composite 
material  being  considered. 

It  is  assumed  that  a relationship  of  the  form  (3)  holds  relating 
the  diffusion  coefficient  (D)  of  the  composite  being  studied  to 
the  temperature  (T) , so  then  for  some  constants  a and  - 


rtct  "Strang  and  G.  Fix,  An  Analysis  of  the  Finite  Element  Method, 
Prentice  Hall,  Englewood  Cliffs,  1973. 


t 

; 


A- 


3 


A 


NSWC/WOL/TR  77-13 


In (D)  = a + B i (A6) 


and  thus  ln(D)  is  a linear  function  of  1/T.  Tne  constants  a and 
B can  be  determined  from  experimental  measurements  such  as  given 
in  Figure  2.  For  example,  for  Narmco  5208/T300' with  a fiber  volume 
fraction  of  .7;  at  T = 348.15  °K  the  value  of  D was  found  to  be 
1.11  E-9  cm2/sec  while  at  T = 303.15°K  the  measured  value  of  D was 
1.21  E-10  cmvsec.  From  these  values  one  determines  that  a = -5198.1 
and  b = -5.6883.  It  is  assumed  that  from  the  weather  data  (or  from 
a kinetic  average  of  the  weather  data)  a subroutine  is  available 
to  obtain  the  temperature  value  (Tn+1)  and  the  relative  humidity 
value  (RHn+l)  to  be  used  at  time  tn+1.  Then 


i 


The  equilibrium  concentration  (g)  of  moisture  at  a given 
relative  humidity  value  (RH)  was  found  by  linearly  interpolating 
experimental  measurements.  The  solubility  of  moisture  in  the  compos- 
ite was  taken  to  be  independent  of  the  temperature  (c.f.  Figure  1). 
For  example,  for  Narmco  5208/T300  with  a fiber  volume  fraction  of 
.7,  linear  interpolation  of  the  following  values  was  used:  (RH=0, 
M=0)  , (RH=80% , M=l.  03%  by^ight  H2O)  , (RH=100%,  M=1.37%)  (c.f. 

Figure  1) . In  this  way  Mq  was  obtained  from  RHn+l.  The  computer 
program  used  to  evaluate  A and  F in  (A5)  is  a much  more  general 
finite  element  code  and  so  is  not  reproduced  here,  as  the  evaluation 
of  the  A and  F given  in  (A4)  is  relatively  straightforward.  An 
example  of  the  determination  of  Tn+1  and  RHn+^  from  given  weather 
data  will  be  discussed  later.  Having  A and  F,  it  remains  to  solve 
the  linear  system  (A5)  to  obtain  Mn+*. 

d . Solution  of  the  Linear  System  AMn+1  = F. 

The  structure  of  the  matrix  A in  equation  (A5)  is  tr idiagonal 
in  that  for  each  row  i of  A,  the  only  nonzero  entries  are  the  di- 
agonal element  Aifi  together  with  Ai_ifi  (when  i/1)  and  Aifi+i  (when 
i/I) . To  solve  AMn+l  = F we  factored  the  matrix  A (by  Gaussian 
elimination)  into  the  form  LU  where  L and  U are  tridiagonal,  L is 
lower  triangular  (i.e.,  Li,j  = 0 when  i<j)  , U is  upper  triangular 
(i.e.,  Ui,j  = 0 when  i > j ) , and  all  the  diagonal  elements  (Lifi)  of 
L are  equal  1 (see  for  example  Section  3.2  of  Chapter  2 of  C5T  ) . 

The  equation  AMn+l  = F is  then  LUMn+l  = F.  Defining  Z = UMn+l,  +. 
one  solves  LZ  = F for  Z,  and  then  one  solves  UMn+l=Z  to  obtain  Mn 
Since  A is  tridiagonal  and  A is  strictly  diagonally  dominant 
(i.e.,  Aj^j  > |a±  j|  for  i=l,...,l),  the  determination  of  L and 

U can  be  done  as  in  the  subroutine  TRIDEC  given  below^.  The 

tsTET  Isaacson  and  H.  B.  Keller,  Analysis  of  Numerical  Methods, 

John  Wiley  & Sons,  New  York,  1966. 


D^tn+1^=  exp  ^ot  + B/Tn+1j, 
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solution  of  LZ  = F and  of  UM  = z is  done  in  the  routine  TRISOL 
given  below.  In  these  routines:  NROWS  denotes  the  number  of  unknowns 
(I);  the  vectors  AL(i),  AD  ( i ) , and  AR(i)  (i  = l, I)  are  defined 
in  terms  of  the  matrix  A by  AL  ( i ) =Aj.-i , i ; AD(i)=Aj  and  AR  ( i ) 

= Ai,i+1  (note  AL(1)  and  AR(I)  are  never  used  in  T&i6ec  or  in  TRISOL); 
RTS I DE  is  the  right-hand  side  (F) ; and  SOLN  is  the  answer  (Mn+1) . 

From  equation  (A4)  it  is  seen  that  the  matrix  A and  vector 
F in  (A5)  depend  on  the  time,  and  thus  we  now  rewrite  (A5)  as 


. n+1 
A 


Mn+1 


,n+l 


(A7) 


Note  that  the  time  dependence  of  An  . is  due  entirely  to  the  D(tn  ) 
term.  Thus  if  D (tn+l) = D(tn)then  An+1  = An  and  so  the  decomposition 
of  An+1  into  the  form  LU  (as  done  in  TRIDEC)  is  the  same  as  that 
for  An  and  hence  need  not  be  repeated.  Thus  the  use  of  a constant 
average  kinetic  temperature  (which  gives  a corresponding  average 
kinetic  diffusion  coefficient)  over  a period  of  time  will  save 
computing  time  because  the  matrix  A need  not  be  recalculated  and 
decomposed  as  long  as  D(t)  stays  fixed  (only  the  formation  of  Fn+1, 
and  the  solution  for  Mn+1  performed  in  TRISOL  need  be  done  at  each 
step) . 


e.  Outline  of  Computer  Program. 

We  now  give  an  outline  of  the  structure  of  the  computer  program 
that  was  used  to  solve  the  diffusion  problem  (7).  The  choice  of 
At  and  the  spatial  mesh  (Aj) , and  an  example  of  a routine  which 
provides  weather  data  will  be  given  below.  The  user  must  supply: 

(1)  a weather  data  subroutine  which  will  produce  the  sequence  of 
temperature  and  relative  humidity  values  to  be  used  - these  could 
be  given  directly  by  meteorological  data  or  by  a kinetic  average 
of  weather  data  (as  discussed  earlier  in  this  report) ; (2)  a diffu- 
sion coefficient  subroutine  which  will  obtain  the  diffusion  coeffi- 
cient of  the  composite  under  study  corresponding  to  any  given 
temperature  value;  and  (3)  a moisture  solubility  subroutine  which 
will  obtain  the  equilibrium  moisture  concentration  corresponding 
to  a given  relative  humidity  value.  The  computer  program  has  a 
SETUP  phase  which  reads  in  necessary  parameters,  and  a MAIN  phase 
which,  given  approximate  solution  values  M.at  time  tn  = nAt, 
produces  the  approximate  solution  values  M7  at  the  next  time  level 
tn+l  = (n+1) At . 1 

Outline  of  Computer  Program 

SETUP 

read  in:  h.  At,  I,  A.  (i=l,...,I),  N = number  of  time  steps  to 
be  taken,  and  interval  of  printout 
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2.  set  Xq  = 0,  Xj  = h,  xi  = x^  + A^  ( i=l, . . . , 1-1) 

3.  read  in  initial  concentration  profile  M°  (i=0,...,I),  m9 
is  the  percent  by  weight  of  moisture  at  point 

4.  check  that  A^  = h 

5.  print  out  the  variables  that  have  been  read  in 

MAIN  at  each  time  step: 

1.  set  tn+^  = (n+l)Afc  and  find  the  temperature  Tn+^  and  the  relative 
humidity  RHn+^  using  the  weather  data  subroutine 

2.  use  the  diffusion  coefficient  subroutine  to  find  the  diffusion 

coefficient  D(tn+^)  corresponding  to  Tn+1  and  use  the  moisture 
solubility  subroutine  to  find  the  equilibrium  moisture  concentra- 
tion corresponding  to  RHn+‘*‘ 

3.  form  the  matrix  An+^  (see  equations  (A7)  and  (A4) ) and  then 

decompose  An+'*'  into  An+"*"  = LU  by  calling  TRIDEC  (skip  this  step 
if:  D(tn+^)  = D(tn)  and  n>0) 

4.  form  Fn+1  (see  equations  (A7)  and  (A4) ) 

5.  solve  An+1  Mn+1  = Fn+1  for  Mn+1  by  calling  TRISOL 

6.  check  if  should  do  printout  of  solution  values  at  tn+1,  if  so, 
printout  tn+"'',  m9+^,  calculate  average  percent  by  weight  moisture 
(10) , etc. 

7.  if  n+l=N  stop,  otherwise  set  n=n+l  and  go  back  to  step  (1)  of 
MAIN 

Samples  of  Computer  outputs  are  shown  in  Figures  18  and  19. 
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